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I- Introduction 

We present preliminary results from a detailed analysis of the X-ray point so urces in 
the XMM-Newton survey of M31 (e.g. iBarnard et aljEobll IPietsch et al.ll2005|) . These 
sources are expected to be mostly X-ray binaries. 

2. Source Detection 

Th e left hand panel of Fig.[l]shows a mosaic of GALEX FUV images of M31 l)Thilker et all 
120051) overlaid by the point sources detected by edetect_chain in SAS 6.0.0. Only sources 
p^j ' which were on both the PN and MOS chips were analysed. The minimum detection like- 
lihood was set to 10 and detections that overlapped at a radius of 20" were excluded. 
Lightcurves, energy spectra and power density spectra were then extracted for each re- 
gion and analysed to determine source properties. Energy spectra were fit with power 
law, bremsstrahlung, black body and neutron star atmosphere models in xspec 11.3.1 
, with free parameters to determine the flux of the source. T he unabsorbed lumino sities 
Q\ ' were then calculated assuming a distance to M31 of 760 kpc l|van den Bergbll2000lL 

^ ■ 3. Results 

In the right hand panel of Fig. ^ we present the first cumulative luminosity functions of 
the XMM-Newton M31 fields to be created using individual fits for each source spectrum 
to obtain the unabsorbed luminosity for each source. The North 1 (solid) and 2 (dashed) 
field luminosity functions are very similar in shape despite having very different numbers 
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of sources. The Core and South 1 luminosity functions have different shapes because these 
fields have not yet been fully surveyed down to the same luminosity limits as the North 
fields. The final luminosity limit for each field will also depend on the good-time interval 
of each observation. We will investigate the significance of these luminosity functions 
■ *h . once the survey is complete. 

We have also created a luminosity function of the North 1 and 2 fie lds combined in or- 
der to compare it to the luminosity function of these fields reported bv lTrudolvubov et al.l 
Unlike in our study, these authors used individual spectral fits only for the bright- 
est sources (Lx > 5 x 10 36 erg s _1 ); the faint sources had their luminosities estimated 
assuming an absorbed power law model with fixed power law index. Above Lx = 10 36 erg 
s _1 our luminosity function can be approximated b y a power law N(L > Lx ) ~ L x a 
with slope a=0.87, inconsistent with a=1.3 found bv iTrudolvubov et alJ ((2002) . This is 
p robably a co n seque nce of the fixed spectral model assumed by these authors. 

iKong et alJ l|2002h used Chandra to study the core region of M31. They extracted 
energy spectra from the longest (8.8ks) observation and carried out spectral fitting only 
on the brightest 20 of the 204 sources they found. All sources were well fit by both 
absorbed power law and Raymond-Smith models. For the other 184 sources, count rates 
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Figure 1. The left hand panel shows the X-ray sources detected, superimposed on a mosaic 
of GALEX FUV observations of M31. The right hand panel shows the luminosity functions of 
each of the fields we have studied. 

were converted to luminosities assuming an absorbed power law model. They obtained 
a luminosity function described by a=1.44 above a break at 1.77 xlO 37 erg s _1 for 
the whole core exposure. This is comparable with our cv=1.96 for 37.7 <log(Lx erg 
s -1 ) < 38.2 in the Cor e luminosity function. 

iGrimm et alJ l)2002j) studied the Milky Way and created luminosity functions for the 
high and low mass binary populations separately. Their cumulative luminosity functions 
had a=0.64 for HMXBs and a=0.26 for LMXBs. Neither of these slopes is consistent with 
our luminosity functions of M31 but the binary population cannot be easily separated in 
M31 and both types of binary are present. 

4. Future Work 

We have so far studied 225 sources over 4 fields and our survey is continuing by 
looking at the unstudied sources we have detected in the 2 outermost disc fields, the 
130 fainter core sources and the 51 fainter sources in the South 1 field. X-ray binaries 
will be identified by their energy spectrum and power density spectrum as outlined in 
iBarnard et al.l l|2004l Eo05|) . We will also continue to create luminosity functions for each 
field to allow us to compare populations. All the luminosity fun ctions we create will 
also need to be corr ected for foreground and bac kground sources llGiacconi et al]l200lt 
lHasinger et al.ll200l[) and for incompleteness fsee lKim fc Fabbianol l|2004l ) 
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